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Abstract 

We derive a "master formula" for the contribution of the three-gluon correlation function 
in the nucleon to the twist-3 single-spin-dependent cross section for semi-inclusive deep- 
inelastic scattering, ep^ — > eDX. This is an extension of the similar formula known for 
the so-called soft-gluon-pole contribution induced by the quark-gluon correlation function 
in a variety of processes. Our master formula reduces the relevant interfering partonic 
subprocess with the participation of the three gluons to the Born cross sections for the j*g — > 
cc scattering, which reveals the new structure behind the twist-3 single spin asymmetry and 
simplifies the actual calculation greatly. A possible extension to higher order corrections is 
also discussed. 
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1 Introduction 



Single (transverse) spin asymmetries (SSA) in high-energy inclusive processes, such as semi- 
inclusive deep inelastic scattering (SIDIS), Drell-Yan process, and hadron production in pp 
collision, appear as a consequence of multi-parton correlations in the hadrons, which did 
not show up in the naive parton picture for high-energy scattering. SSA reveals the multi- 
parton correlation most unambiguously as a leading twist-3 observable in the framework of 
collinear factorization which is valid to describe hadron production with large transverse 
momentum ^> Aqcd- 1 Among those correlations, the effect of quark-gluon correlations 
have been widely studied in the literature and our understanding on the mechanism of 
SSA has made a great progress [1]-[16]. There are also some studies on the contribution to 
SSA from the purely gluonic correlation in the nucleon [17, 18, 19, 20] and the multi-parton 
correlations in the fragmentation functions [21, 22]. 

In our recent paper [19], we have established the formalism for calculating the twist-3 
single-spin dependent cross section induced by the three-gluon correlation functions in the 
transversely polarized nucleon. There, we clarified a complete set of the gauge-invariant 
three-gluon correlation functions and derived the corresponding cross section for the IP- 
meson production in SIDIS, ep^ — > eDX, which is relevant to probe the purely gluonic 
correlation leading to SSA. 2 In particular, we have proved the factorization property of 
the twist-3 cross section and given a detailed prescription of expressing the cross section in 
terms of the gauge-invariant three-gluon correlation functions. Our result differed from the 
previous study [17, 18], and we clarified why the previous result needs to be corrected. 

The partonic cross section in this formalism is given as a pole contribution of an internal 
propagator in the hard part, reflecting the naively T-odd nature of SSA. The pole forces 
one of the gluon lines in the three-gluon correlation functions to be soft, and hence the 
pole is called the soft-gluon-pole (SGP). For the SGP contribution from the quark-gluon 
correlation functions, it has been shown in [10, 11] that the corresponding twist-3 hard cross 
sections for ep^ — » enX and p^p — > ttX have a simple relation with the twist-2 unpolarized 
hard cross section for the same processes. This connection greatly facilitates the actual 
calculation and makes the structure of the SGP contribution transparent. 

The purpose of this paper is to show that the whole three-gluon contribution to ep* — > 
eDX can also be obtained from the Born cross sections for the 7*g — » cc scattering. In 
particular, we will show that some of the twist-3 cross section is completely determined by 
the gluon-contribution to the twist-2 unpolarized cross section. 

The remainder of this paper is organized as follows: In section 2, we briefly recall a 
complete set of the twist-3 three-gluon correlation functions in the transversely polarized 

1 For the contribution from the quark-gluon correlation functions, consistency between this approach 
and the transverse-momentum-dependent factorization in the region Aqcd <C Pt *C Q with the hard scale 
Q has also been shown for some processes [6, 7, 8]. 

2 Thc D-meson production at large Pt in SIDIS is dominated by the photon-gluon fusion subprosess. 
Other possible sources of the asymmetry in the D-meson production are associated with the charm quark 
content in the polarized nucleon, such as the c-quark transversity distribution and the c-quark-gluon corre- 
lation functions, but those intrinsic charm contributions are expected to give tiny corrections in the large 
Pt region which we are interested in this paper. 
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nucleoli. In section 3, we summarize the twist-3 formalism for the three-gluon correlation 
functions developed in [19]. We mostly follow the notation used in [19], but introduce a 
slightly different convention for the azimuthal angles and the hard part, which turned out 
to be more convenient for our purpose. In section 4, we develop a new master formula which 
connects the three-gluon contribution to the twist-3 single-spin-dependent cross section to 
a simpler cross section for the 7*g — > cc scattering. We also discuss the extension of the 
formula for higher order corrections. Section 5 is devoted to a brief summary of the outcome 
of this paper. 

2 Three-gluon correlation functions in the transversely- 
polarized nucleon 

As clarified in [23, 24, 19], there are two- independent twist-3 three-gluon correlation func- 
tions in the transversely-polarized nucleon, 0(xi,X2) and N(xi,x 2 ), which are the Lorentz- 
scalar functions of the two momentum fractions x\ and x 2 , defined as 

O a ^(x u x 2 ) = -gi 3 f—f ^e^e^-^ipSld^F^^F^^F^iXn^pS) 
J 2ix J 2iv 

= 2iM N [0{x l ,x 2 )g a ^e^ nS + 0{x 2 ,x 2 - x^g^e^ + 0(x 1 ,x 1 - x 2 )g^ a e ppnS ] , (1) 

N a ^( Xl ,x 2 ) = -gi z f—f ^e tX ^e i ^- x '\ P S\if bm F^ n ^)F2 n {iin)F^\\n)\ V S) 
J 2tc J 2ir 

= 2iM N [N^x^g^e 1 ^ - N{x 2 , x 2 - x 1 )g^e apnS - N{x 1 ,x 1 - x 2 )g^ a e^ nS ] , (2) 

up to the irrelevant terms of twist-4 and higher, where F" n = F^rip with F^ = d a A^ — 
<9 /3 A" + gf abc A^A^ being the gluon field strength tensor, d bca and f bca are, respectively, the 
symmetric and anti-symmetric structure constants of the color SU(3) group, and we have 
suppressed the gauge-link operators which appropriately connect the field strength tensors 
so as to ensure the gauge invariance. p is the nucleon momentum, S is the transverse spin 
vector of the nucleon normalized as S 2 = — 1, and Mjv is the nucleon mass so that 0(x±, x 2 ) 
and N(xi,x 2 ) are dimensionless. In the twist-3 accuracy, p can be regarded as lightlike 
(p 2 = 0) and n is another lightlike vector satisfying p ■ n = 1. To to be specific, we take 
pt 1 = (p + ,0,0j_), n M = (0,n~,0_i_) and = (0,0, SjJ. Hermiticity, invariance under the 
transformations P and T, and the permutation symmetry among the participating three 
gluon-fields imply that 0(xi,x 2 ) and N(xi,x 2 ) are real functions and satisfy the relations, 

0(X!, x 2 ) = 0(x 2 , Xi), 0(X!, x 2 ) = 0(-x 1 , -x 2 ), 

N(x 1 ,x 2 ) = N(x 2 ,x 1 ), N(x u x 2 ) = -N(-x 1 ,-x 2 ). (3) 
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3 Summary of the twist-3 formalism for ep* — » eDX 
3.1 Kinematics 

Here, we summarize the kinematics for the SIDIS process, 

e(£) + pt(p, S) e(f ) + + X, (4) 

where the final D-meson has the mass rrih, i.e., = m^. This process is described by the 
five independent Lorentz invariants: 

S ep = (p + £)\ x b3 = Q 2 = -g 2 = -{I - £') 2 , */ = q T = f^!, (5) 

2p ■ q P ■ Q 

where q t is the "transverse" component of q defined as 

qt ~ q + {(P-P h ) 2 P-P h ) P P-P h h ' ( ) 

satisfying q t • p = q t ■ P h = 0. In the actual calculation we work in the hadron frame where 
the virtual photon and the initial nucleon are collinear, i.e., both move along the z-axis. In 
this frame, their momenta q and p are given by 

«" = (?°,9) = (0,0,0,-Q), p" = ^,0,0,^ . (7) 

The azimuthal angle of the hadron plane as measured from the xz plane is taken to be x 
and thus the momentum of the D-meson is parameterized as 

p M ZfQ m h 2( 1t 2g T . q\ m\ \ 

Ph = — [ 1 + ^ + Z jQ>^ COSX ^ SmX >- 1 + ^ + ^) ■ (8) 

The transverse momentum of the D-meson in this frame is given by PhT = Zfqr, which is 
true in any frame where the 3-momenta q and p are collinear. The azimuthal angle of the 
lepton plane measured from the xz plane is taken to be and thus the lepton momentum 
can be parameterized as 

£ M = — (cosh^, sinh ip cos 0, sinh -0 sin 0,-1) , 

= — (cosh -0, sinh -0 cos 0, sinh -0 sin 0, 1) , (9) 
2 



where 



cosh^ = ^^-l. (10) 

Q 
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We parameterize the transverse spin vector of the initial nucleon S p as 



^ = (0, cos $ 5 , sin $ 5 ,0), (11) 

— * 

with the azimuthal angle $5 of S. Although three azimuthal angles 0, x and $s are defined 
above, it is obvious that the cross section for ep^ — > eDX depends on them through only 
the relative angles — x an d $s — X- Thus, it can be expressed in terms of S ep , Xbj, Q 2 , 
Zf, gf., — x an d $s - x in the above hadron frame. Note that 0, x an d &s are invariant 
under boosts in the ^-direction, so that the cross section presented below is the same in 
any frame where q and p are collinear. 

With the kinematical variables defined above, the differential cross section for the D- 
meson production in SIDIS using the unpolarized lepton can be calculated with 



d 6 a a 2 em 



dxbjdQ 2 dzfdq^d(pdx 1287r 4 a;^S'| p (5 



z f Lr{l,e)W^[p,q,P h ), (12) 



where L^ u (£, £') = 2(£ p £' u + £ u £' p ) — Q 2 g pu is the corresponding leptonic tensor, W^ u (p, q, P h ) 
is the hadronic tensor in the same normalization as in [19], and a em = e 2 /(4-7r) is the QED 
coupling constant. The single-spin (S) dependent part in the cross section (12) describes 
the SSA in the process (4). For this part, one can transform the azimuthal element in the 
LHS of (12) as d<pdx — > d(j)d<& s and set x = 0. In fact, in our previous paper [19], and 
$5 were used as the azimuthal angles, respectively, for the lepton plane and the spin vector 
measured from the hadron plane by setting x = from the beginning. In the expression 
for the differential cross section derived in [19], the differential element "gM>s" appearing in 
the above sense was missing, and thus should be supplied together with the factor 1 / (2n) 
for the cross section with all the results unchanged. 



3.2 Twist-3 cross section 

As shown in [19], the contribution from the three-gluon correlation functions to the hadronic 
tensor is relevant for having the sizeable SSA in the collinear factorization to describe 
the large- Put -D-meson production and is represented by the diagrams of the type shown in 
Fig. 1. Here, the twist-2 fragmentation function D(z) for a c-quark to become the D-meson 
is already factorized as the upper blob. In the other part of the diagrams, the partonic hard 
part S^ / c pTX (ki,k 2 ,q,p c ), represented by the middle blob, is combined with the correlation 
functions ~ (AlA^A p a ) of the gluon fields, A p a {£), in the nucleon (lower blob), where p, r, A 
and a, b, c are, respectively, Lorentz and color indices for the relevant three gluon-fields. 
fi and v in S°i^ c pTX (ki,k2,q,Pc) represent the Lorentz indices for the virtual photon. Com- 
pared to the conventional photon-gluon fusion subprocess, 7*5 — > cc, an additional gluon 
participates in the partonic hard part S^^^ki, k 2 , q,p c ) and allows us to obtain interfering 
phase arising from the unpinched pole contribution of an internal propagator contained in 
x (ki,k2,q,p c ). Owing to the symmetry property of the spin-dependent part of the 
nucleon matrix elements {AlA^A p ) under the P- and T-transformations [19], only such 
interfering contribution in the hard part S^: pTX (ki, /c 2 , q,p c ) can give rise to the single-spin- 
dependent cross section. In the leading order with respect to the QCD coupling constant, 
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the corresponding contribution to x (k±, k 2 , q,p c ) can be obtained from the diagrams 

shown in Fig. 2 (together with their mirror diagrams), where the bar on a propagator 
indicates that the pole part is to be taken from the propagator; in principle, other prop- 
agators can produce similar pole contributions, but we need not consider separately those 
contributions, which indeed correspond to the diagrams obtained by the permutation of the 
gluons in the diagrams of Fig. 2. It has been shown in [19] that the total contribution to the 
single-spin-dependent cross section from Figs. 1, 2 can be expressed in a gauge-invariant 
form as 



W^(p,q,P h ) 



dz 



D(z) 



dxi 

Xi 



dx 2 

X 2 



X 



dS^ pTX {k 1 ,k 2 ,q,p c )p 
dk% 



pole 



up to the twist-3 accuracy in the collinear factorization with k\ 

P — „P _r,P^ \A a ^ 



u p a uj T f3 uj^M a F % c (x 1 ,x 2 ), (13) 
x\p and k 2 = x 2 p, where 



U) 



9 r a —p^na and M.p] bc {xi, x 2 ) denotes the three-gluon lightcone correlation functions 



defined in terms of the gluon field-strength tensors, as 

d\ f dfj, 



M a F % c (x u x 2 ) = -gi 3 J^j ^e^ 1 e^--)(p5|Ff"(0)F^(HFr(Ar i )|p5) 



with O a ^{x\,x 2 ) and N al3l (xi, x 2 ) given in (1) and (2). We introduced the notation 
where 



(14) 

ipolc 



implies that the bared propagator arising in Fig. 2 (and 



d be replaced by its pole part oc 5 {{p c + k x — k 2 )' 



or 



S Zi c ,pT\( k i,k 2 ,q,Pc) 
its mirror diagrams) shou 
5 ((p c + k 2 — ki) 2 — ml)), with p c the momentum of the c-quark fragmenting into the D- 
meson, p\ = m 2 ,; this fixes a momentum fraction at x\ = x 2 in the collinear limit ki = xip 
in (13). p c is parameterized by the momentum fraction z associated with the fragmentation 
function D(z), as 



z h ^2p-P h 



m 2 c z 



mt 



2 I Q 2 z^Q 2 ' 



2q T 2q T . 

- — cosx, -p-smx, 



Q 



Q 



2 2 ' 

+ q 2 + 2 2 g 2 



(15) 



where z = Zf/z, and the second line shows the explicit form in the hadron frame with (7), 
(8). The summation over the c and c quark contributions, as well as the corresponding flavor 
index on D(z) and x (ki, k 2 , q,p c ), is implicit in (13). Actually, many terms of twist-3, 
other than those in (13), are generated by the collinear expansion of the diagram of Fig. 1, 
and appear to be gauge- noninvariant. We emphasize that all these gauge-noninvariant 
twist-3 terms either cancel or vanish in the cross section, as demonstrated in [19] by the 
use of Ward identities. 
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Figure 1: Generic diagram giving rise to the twist-3 contribution to the hadronic tensor of 
ep^ eDX induced by the gluonic effect in the nucleon. It is decomposed into the nucleon 
matrix element (lower blob), .D-meson matrix element (upper blob), and the partonic hard 
scattering part by the virtual photon (middle blob). 

3.3 Calculation of L^W^ 

To calculate the contraction L^ u (£, i')W^ v (p, q, Ph) in (12), we introduce the following four 
vectors which are orthogonal to each other, similarly as in [19]: 

= e» vpc Z v X p T a , 

* = (16) 

In the hadron frame specified by (7) and (8), these vectors become = (1,0,0,0), X^ = 
(0, cosx, sinx, 0), F M = (0, — sinx, cosx, 0), = (0, 0, 0, 1), so that 

(0,1,0,0) =cos X ^-sin X ^, 

(0,0,1,0) =sin X ^ + cos X ^, (17) 

In the present case, W^ u can be expanded in terms of the following six independent ten- 
sors [19]: 



7 




Figure 2: Feynman diagrams for the partonic hard part in Fig. 1, representing the photon- 
gluon fusion subprocesses that give rise to the pole contribution for ep^ — > eDX in the 
leading order with respect to the QCD coupling constant. The short bar on the internal 
c-quark line indicates that the pole part is to be taken from that propagator. In the text, 
momenta are assigned as shown in the upper-left diagram, where p c denotes the momentum 
of the c-quark fragmenting into the Z)-meson in the final state. The mirror diagrams also 
contribute. 

vr = t»x u + x^t u , vr = x ^ xu - y ^ y ^ 

Vg v = T^Y U + Y^T U , Vr = X"Y" + Y»X V . (18) 
We also introduce the inverse tensors V£ v for the above V^ v : 

VT = ^(2T^T u + X^X V + Y^Y U ), V% v = T^T", 

Vf = --(T fi X v + XT"), Vf = -{X^X" - Y^Y U ), 

T>r = —(T»Y V + FT), 9^ = ~{X»Y V + W). (19) 
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With these definitions, one has 



E [L^vn{w P „vr}=Q 2 E A(0-x) 

k=l, -,4,8,9 fc=l, ■■■,4,8,9 



(20) 



where — x) = L^V^ /Q 2 parameterize the dependence of the cross section on the 

azimuthal angle of the lepton plane relative to the hadron plane (see (9), (8)), and depend 
on and x through — x, with 



AM = 


1 + cosh 2 -0, 


A 2 (cf>) = 


-2, 


-4 3 (0) = 


— cos0sinh 2-0, 


M4>) = 


cos20sinh 2 -0, 


M<t>) = 


— sin0sinh2 , 0, 


M4>) = 


sin 20sinh 2, 0. 



(21) 

By the expansion (20), the cross section for ep^ — > eDX consists of the five structure 
functions associated with Aip, A3, A±, As and Ag, respectively, which have different de- 
pendences on the azimuthal angle 0. 



4 Master formula for three-gluon contribution 

4.1 Connection between the twist-3 hard part and the j*9 ~~ ^ °£ 
scattering 

To obtain the twist-3 cross section based on (13), one has to calculate the corresponding 
hard part as the derivative, dS^ c a/3X (ki,k2,q,p c )p x /dk2, for the contributions from the 
diagrams in Fig. 2, followed by the collinear limit ki — > Xip. We note that many building 
blocks in the diagrams in Fig. 2 depend on fc 2 due to the participation of three external 
gluons with the momenta hi, k 2 and k 2 — ki, so that taking the derivative with respect to k\ 
produces many terms in the intermediate step and is quite complicated. As will be shown 
below, the above hard part for the twist-3 cross section is connected to the hard part with 
only the two external gluons, representing the "twist- 2 level" partonic scattering, 7*^ — > cc. 
This implies, in particular, some relevant contribution to the hard part for the twist-3 
cross section is completely determined by the hard part for the twist-2 unpolarized process 
ep — > eDX, which is associated with the gluon density-distribution in the nucleon, and, 
similarly, the spin- dependent contributions in the partonic scattering 7*^ — > cc completely 
determine the remaining hard part for the twist-3 cross section. 

In order to prove this, we proceed similarly as in the proof in [10, 11] for the mas- 
ter formula associated with the quark-gluon correlation, but we discuss each step of our 
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proof in detail because it involves extensions for the case not only with the three-gluon 
correlation, but also with the nonzero quark-mass, compared to the massless quark case 
treated in [10, 11]. We first note that the hard part shown by the diagrams in Fig. 2 has 
the structure obtained by attaching an extra gluon-line to the c-quark line in the diagrams 
in Fig. 3, which show the leading-order contribution for the 7*g — > cc scattering with the 
c-quark fragmenting into the final D-meson. We denote the sum of the contributions of the 
diagrams in Fig. 3 as Sj%}f£p(xp,q,p c ), where the Lorentz indices a and f3, as well as the 
color indices a and b, are associated with the external gluon lines that have the momentum 
X P- S^^pixp-, q,Pc)/8 with the color indices averaged over represents the hard part for the 
twist-2 cross sections. In particular, the twist-2 unpolarized cross section in (12) is given 
by the contribution to the hadronic tensor, 

W£(P, ?, Ph) = I %m f ^S^(xp, q, Pc )G a J(x), (22) 

where Qab( x ) is the light-cone correlation function of the gluon's field strength tensors in 
the nucleon, defined as 

Q ^ {X) = ll ^ iXx ^\F b ^(0)Fr(\n)\pS) = -\gf x \b ah G{x) + • • • , (23) 

with the unpolarized gluon density G(x), and g^ 3 = g alS — p a n^ — n a pP . From Fig. 3, 
^m^a^ViQ.iVc) is obtained as 



S jS?ap( x P> Q, Pc) = Tr [Fp(xp, q, p c ) (tf c + m c ) T^(xp, q, p c )V(xp + q- p c )\ , (24) 

where Tr[- • •] indicates the trace over both Dirac and color indices, {yj> c + m c ) with (15) is 
the projection matrix for the D-meson fragmentation function, and the factor, 

V(k) = 2tv(^ -m c )5(k 2 -m 2 c ), (25) 

is associated with the final-state cut for the unobserved c-quark with the momentum k 
flowing from the left to the right of the cut. The factor J = a{ x P->1^Pc) is the 7gcc-vertex 
function in the left of the cut, containing the photon-quark and the gluon-quark vertices 
linked by the quark propagator; here, the Lorentz and color indices, a and a, are associated 
with the external gluon, while the Lorentz index \x for the virtual photon is suppressed 
for simplicity. The similar factor in the right of the cut can be written as Fp(xp,q,p c ) = 
7 °[^(xp,g,p c )]t 7 o. 

With the replacement xp — > k\ in (24), where k\ has nonzero transverse components, 
s( $?ap( x P, Q,Pc) may be extended to S { ^{k u q,p c ). Sffi^ki, k 2 , q,p c )p X arising in (13) is 
obtained by attaching an additional gluon to S^J a ^(ki,q,p c ), in particular, to the c-quark 
line fragmenting into the D meson, where the gluon carries the momentum ki — hi, the 
polarization associated with p A , and the color index c. When the extra gluon is attached 
in the left of the cut as in Fig. 2, we obtain 



Sl%(ki,k 2 ,q,p c ) =Tr j£(fc 2 , q,p c )t c L(h - k 2 + p c )KiK q, h - k 2 + p c )V(k 2 + q - p c 



(26) 
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.rp 







Figure 3: Leading-order diagrams for the 7*g — >■ cc scattering cross section. 



where 



L(ki - k 2 +p c 



+ m c 



-1 



m c + ie 



(27) 



and t c is the color matrix for the quark-gluon vertex (the coupling constant g for the 
quark-gluon vertex is absorbed in A^^^ c (xi, x 2 ) of (14)). Likewise, when the extra gluon 
is attached in the right of the cut as in the mirror diagrams of Fig. 2, we obtain 



S a R %(kl,k 2 ,q,p c ) = Tr [J* p (k 2 , q, k 2 - k x + p c )t c R(k 2 - h + p e )J*(k 1 ,q,p c )V(k 1 + q - p c ) 

(2f 

with 



R(k 2 - fa + p c 



m r — ie 



+ m c )= 1 L\k 2 -fa+p c ) 1 °, (29) 



so that 



SZp(k u k 2 ,q,p c ) = S b L % (k 2 ,fa,q,p c ) 



Thus, S°£f. a a X (ki, k 2 , q } pc)p can be obtained as the sum of (26) and (28): 



S^pxiK k 2 , q, Pc )p x = S a L %(h, k 2 , q,p c ) + S a R %(h, h, q,p c ). 



(30) 



(31) 
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Comparing (24) with (31), we see that Sff: p X (ki, k 2 , q,p c )p x is obtained from S ( ^J^ /3 (k 2 , q,p c ) 
by the formal replacements, {jf> c + m c ) — > t c L{ki — k 2 +p c ) and (#i c + m c ) — > t c R(k 2 — ki+p c ), 
together with the appropriate momentum-shifts in the remaining factors. 

We calculate 9S^ a/3X (ki, k 2 , q,p c )p x /dk 2 \ ki=x _ p in (13) based on (26)-(31), keeping their 

relevant structure that manifests the above-mentioned correspondence with S ( / ^ a c ^ /3 (k 2 , q,p c )- 
For this purpose, we need the collinear limit of (27) and (29), 



L ((xi - x 2 )p + p c ) = -R ((x 2 - xi)p + p c ) 

and of their derivatives, 
(9L(A;i - k 2 +p c 



+ m c ) 



-1 



%i — x 2 + ie' 



(32) 



dk% 

dR(k 2 - kx + p c ) 



(|( c + m c )^7 a 



+ m c 



dk? 



ki=Xip 



7 v>,(33) 

2p ■ p c xi — x 2 + ie p ■ p c (xi — x 2 + ie) 1 



+ m c ) 



2p-p c 



xi — x 2 + 



+ 



Pea (^ c + m c) 



1 



P ■ Pc (%i ~x 2 + ief 



.(34) 



In these relations (32)-(34), the poles relevant to [■ ■ -] polc in (13) are unveiled; it is straight- 
forward to see that the 7gcc-vertex functions and their derivatives do not produce the pole 
contributions in the ki — > XiP limit with Xi > 0. Based on the pole structures in (32)- (34), 
we evaluate the pole contribution in (13) as 



dsf^jk^q^cy 



dk2 





pole 


ki — %iP_ 





ki=Xip 



(i)+(ii)+(iii) 



(35) 



decomposing it into three parts (i)-(iii), where (i) denotes the pole contributions from the 
second term in (33) and (34), (ii) denotes those from the first term in (33) and (34), and 
(iii) denotes the remaining pole contributions due to (32). We obtain 



dkl 



ki=Xip 



(i) 



Pcy 
P-Pc 



- pole 



(x 1 - x 2 + ie) 2 



xTr 



-T\{x 2 p, q, p c )t c (J c + mc) Fa{xip, q, (x 1 - x 2 )p + p c )V(x 2 p + q - p c ) 



+F b Jx 2 p, q, (x 2 - Xi)p + p c )t c i$ c + m c ) T^xtf, q, p c )V{x x p + q - p c ) 



Per/ 
P-Pc 

xTr 



1 



x\ — x 2 + ie 



pole 



d 



Fp(xiP, q,Pc)t c ($ c + m c ) I p^-^Taixip, q,p c ) j V{x x p + q-p c 
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+ (p^-^J^ix^p, q,p c )^j t c (tf c + m c ) Tlixxp, q,p c )V[x x p + q-p c ) 
+J r p(x 1 p, q,p c )t c (jzf c + m c ) T a a {x x p, q,p c ) (p^-^Vixxp + q - p c ) 



Pc^f 
P-Pc 

xTr 



-1 



pole 



X\ — x 2 + ie 
d 

P^q^^^P^ 9' P^ $c + m c) Ki x iP, 9, P^fD{xYP + q-p c ) 



-Tl{x x p, q, Pc)t c p'J r ^(x 1 p, q, p c )V(x l p + q~p c ) 



(36) 



In the second equality in (36), we have performed the Taylor expansion with respect to x 2 
around x\ for the contributions inside the trace Tr, and the partial derivative with respect 
to p c implies the shorthand notation of 



d 



dp. 



Tlf(Pc) 



d 



f(r) 



(37) 



r-s>p c 



for a function /(r) of a four- vector r M . The leading-order (zeroth-order) term in the Taylor 
expansion would give rise to the double pole contribution to (36), but the corresponding 
contributions generated from the two terms in the RHS of (31) cancel; this type of cancella- 
tion may be considered as a result of gauge invariance [10] . On the other hand, the second 
and higher-order terms in the expansion give the vanishing contribution to (36) when com- 
bined with the factor [l/(x\ — x 2 + ie) 2 ] polc . Similarly, for the contribution labeled by (ii) 
above, we easily obtain 

(") 



dS^ a s X {ki,k 2 ,q,p c )p x 



dkl 



x\ — x 2 + ie 



pole 



Tr 



Pay 



P-Pc 

The contribution labeled by (iii) can be expressed as 

(iii) 



4 ^(xip, q, p c )V(x 1 p + q-p c ) 



.(38) 



dS^ afjx (k u k 2 ,q,p c y 



dkl 



ki=Xip 



xi — x 2 + ie 



pole 



Tr 



d ^ 



dp 



-T^XxP, q, p c )t c i$ c + m c) K( x iP, 9, Pc)1>{xiP + q-p c ) 



-Tl (x!p, q, p c )t c 7 7 T^Xtf, q, p c )V{x x p + q - p c ) 



(39) 
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where we have replaced the derivative d/dk 2 by the relevant d/dpj, and have set x 2 = x\ 
inside the Tr due to the presence of the factor [\/{x 1 - x 2 + ie)] polc . The sum of the results 
(36)- (39) yields the compact form for (35), 



dkl 



pole 



— in5(xi — x 2 ) 



dpi p ■ p c dpc 



where 



(40) 



Sfu-^ixiP.q^Pc) = Tr[^{x 1 p,q,p c )t c ^ c + m c )^{x 1 p,q,p c )V{x 1 p + q-p c )] . (41) 
Noting that the RHS of (40) vanishes when contracted by p y , (40) may be calculated as 



dkl 



pole 



-inSix! - x 2 )-^Sf u c . ap {x l p ) q,p c ) ) (42) 



where the derivative is to be taken under the on-shell condition p c = m c , regarding pj as 
a variable dependent on p~ and pi' 2 as in (15). Comparing (41) and (24), S°^ a p(xip, q,p c ) 

is the same as S^j a ^(xip,q,p c ), except that the former has an extra insertion of the color 
matrix t c . The relations (40) and (42) clearly indicate that the three-gluon contribution to 
the twist-3 cross section for ep^ — > eDX (see (13)) can be derived from the cross sections 
for the 7*g — > cc scattering, and we find that those relations, obtained for the three- 
gluon contribution and with a massive quark, have the structure formally similar to the 
corresponding relations derived in [10, 11] for the quark-gluon contribution and with the 
massless quarks. Substituting (13), (42) into (12), we obtain the master formula, the main 
result of this paper, which allows us to derive the explicit form of the whole contribution 
to the twist-3 SSA in ep t ->■ eDX, as we demonstrate in the next section. 



4.2 Calculation of the twist-3 cross section based on master for- 
mula 

We are now in a position to carry out the derivative arising in (42). In the formula (13) 
with (14) and (42), it is sufficient to consider the corresponding derivative for the color- 
averaged components of Sffi: a p(xp, q,p c ), which coincide with the color- averaged component 

°f ^^a^PiliPc), as apparent from the color structure in (24), (41). We thus define the 

color-averaged hard part S^J. a/3 (xp,q,p c ) as the appropriate color contraction: 

\^S<2*(xp, q, Pc ) = \s^(xp, q,p c ). (43) 

Then, we have 

^d abc S^(xp, q,p c ) = ^f abc S; b ^(xp, q,p c ) = \s^(xp, q,p c ). (44) 
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Furthermore, the above derivative can be calculated most conveniently after taking the 
contraction of S^J. a/3 (xp,q,p c ) with the leptonic tensor L^, as implied in (12). According 
to (20), we introduce 

VFsW^xp, q, p c ) = 2ir—H*p(xp, q, p c ). (45) 

Here, in the RHS, we separated the coupling constant, contained in the jgcc- vertex function 
FaixPiQiPc) in (24), as well as the factor 2tc/(zQ 2 ) contained in V(xp + q — p c ) for the 
unobserved final-state, such that (see (25)) 



2/02 



^ 2 g 



(46) 



with x = xtj/x and £ = Zf/z. Thus, one obtains 



ITW^ = 2 -^\ £ / ~D(z) I %{-™)^ [A k (<t>- X )H k a ,(xp,q,p c 

z f 4 fc=l, --,4,8,9 J Z J X a P c 

x (o a ^(x, x) + A^x, x)) , (47) 
for the single-spin-dependent contribution, where (see (1) and (2)) 

O a ^(x,x) = 2%M N [0(x,x)gfe^ pnS + O(x,0)(gi J e apnS + gTe^ pnS )] , 

N a ^(x, x) = 2iM N [N(x, x)gfe 1pnS - N(x, 0){g^e apnS + gJ"e^ nS )\ , (48) 

up to the irrelevant terms corresponding to the twist higher than three. We note that (47) 
with (48) is described by the two kinds of partonic hard parts: the hard part associated 
with 0(x,x) is same as that for N(x,x), but 0(x, 0) as well as N(x,0) accompanies the 
hard part of another kind. 

The result (47) shows that only the derivative with respect to the transverse compo- 
nents, (Pc)Pc) = Pc-L) contributes to the twist-3 cross section. As noted below (42), those 
components can be varied as independent variables when performing the derivative, and, 
based on the representation in (15), the corresponding derivative may be performed through 
the magnitude p c ± = |p c j_| = zqx and the azimuthal angle x °f the transverse components, 



as 

d d sinx d d . d cos x d , . 

7 - 1 = cosxj: ^l = sm *x + 7T- ( 49 ) 

OPc °Pe± Pel. OX OPc °Pc± Pel. OX 

The formulae (21) indicate that d/dp c ± hits only H^(xp, q,p c ) in (47). Substituting 
(48) into (47) and using e lpnS = sin $5 and e 2pnS = — cos $5, the hard cross section for 
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{0(x, x), N(x, x)} can be obtained as 



dp, 



c 



= Sing^-x) _ g^gf + co s( . s - x) dM 4, - x) f 

where we have used the fact that the scalar function H^ag^ does not depend on the angular 
variable x- Indeed, comparing with (22), (23), we remark that —H^g°^ is nothing but 
the twist-2 partonic part for the unpolarized gluon density G(x), and its explicit form is 
calculated in [19]. Accordingly, the partonic hard part for {0(x, x), N(x, x)} is completely 
determined from the knowledge on the twist-2 unpolarized cross section. Using the rela- 
tions = —A$ and = — 2.A9, we find that the terms associated with the azimuthal 
structures .4.8,9 arise from the second term of (50), accompanying cos(<E>s — x), although 
such azimuthal structures were absent from the twist-2 unpolarized cross section (22), i.e., 
H^g°f = for k = 8,9 (see (56) below). These remarkable features revealed in (50) have 
been observed similarly in the master formula for the SGP contribution to the SSA induced 
by the quark-gluon correlation [10]. We also remind that the result (50) depends on the 
azimuthal angles through <fi — x and $5 — Xi as noted in section 3.1. 3 

The hard part for {0(x, 0), N(x, 0)} can be calculated similarly. In this case, however, 
we encounter the non-scalar components of H*p, which have the dependence on the angular 
variable x- The x dependence of this type was irrelevant for the cases discussed in [10, 11], 
but gives the novel contribution in the master-formula approach for the present case. In 
order to deal with those non-scalar components, we write down the corresponding hard 
part explicitly as 



A { A(0 - x)H^} (g^ nS + gTe^ s ) 



-2^4^ - x)H k u e lpnS - 4xM4> - X){H k l2 + H* x )f*»* 
dpi dp\ 

9 r A(0 - x)(H k 12 + H k 2l )e^ nS - 2^A k (<f> - X )H k 2 e^ nS 



= -2 sin $ 5 fcos X /- - Ak (<P - x)H k n 

V dp c± p cL dx 



3 Since this dependence is obvious from the beginning, one may alternatively consider the derivative in 
(47) with (49) at \ — to get ^-dependence and restore the x-dependence by the shift $5 — > $5 — \ 
to reach (50), which is much simpler. In other words, one can calculate the cross section for \ ~> by 
regarding \ as the small parameter that allows us to take the derivative d/dp2 in (47). This also applies 
to (52) below. 
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+ COS $5 cos x 



d 



dp c ± Pel. dx 



21) 



- sin $ 5 (sin X -A- + C ^A®\ M4> - x)(H k 12 + H k 2l ) 
V dp c± p c± dx) 

( ■ 9 cos x d \ . , , . TTk 
+2 cos $ 5 sin Xt; + -jr- A k - X # 2 2 , 



(51) 



and thus we need the x dependence of H^, H\ 2 + if 21 and H 22 - This dependence can be 
separated using (17) like H\ 2 = H*p (cos xX a — sin x^ a ) (si 11 + cos xX^) '■> as apparent 
from the definition (16) for the basis vectors constructed with the momenta that are asso- 
ciated with the hadron plane, H xx = H^X^X", H XY = H^X^Y", etc, are unchanged by 
the rotation of the hadron plane around the z-axis, i.e., are independent of x- After some 
algebra, one obtains, in the obvious notation, 



d 
dp, 



{ A(0 - x)H k aP ) (g^e apnS + gT^ 8 ) 
■2A k ( 



= sin($ 5 - x) 



- X) dH k xx 

z dqx 

1 ( dA k {ct> - x) uk 



zq T 



H(XY+YX) + 2-4fe(0 — X)H(_XX+YY) 



+ COs($5 - X) 



A(0 - X) dH (XY+YX) 

z dqT 

dA k ((j)-x) rjk 



zq T 



-2 



d<j) 



H* Y + 2A k (<f>-x)Hf XY+YX ^ 



(52) 



To proceed further, we define the partonic hard cross sections for the "7*^ — > cc" scattering, 
separating the delta function of (46) and the electric charge e c = 2/3 of the c-quark from 
the hard parts arising in (50), (52): 



-H k pgf = elaV(Q,q T ,x,z)6^ 2 - 
-H k xx = ela^(Q,q T ,x,z)5^- 



-H YY = ejal Y (Q,q T ,x,z)5 



,Q 2 



1 - - ) ( 1 " t ) + 

x 



1 - 



X 



X 



1 - 



1 - 



+ 



+ 



2D2 



z 2 Q 
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-H, 



(XY+YX) — e c a k 



2~{XY} 



(Q,q T ,x,z)S 



,Q 2 



+ 



m c 



(53) 



where the partonic cross sections a J k (j = U, XX, YY, {XY}) are the functions of Q, qr, x 
and z. When (50) and (52) using these forms are inserted into (47), the derivative d/dqr 
hitting the delta function in (53) can be treated by integration by parts with respect to x; 
for example, the contribution from the first term in the RHS of (50), convoluted with the 
correlation function 0(x,x), can be calculated as 



1 d 



zdqx 
2qr 

" Q 2 



dx 



°k(Q,<lT,x,z)0( y x,x)5 



x< 



dx 
x 2 



X&k(Q,<lT,X,z) 



1-z 



x- 



Q 2 

dO(x, x) 



- 1- 



1 



X 



1 - 



1 



+ 



to; 



dx 



20(x,x) 



z 2 Q 2 J 
Q 2 da%(Q,q T ,x,z) 



dq 



x 2 da^(Q,q T ,x,z) 



1-z 



dx 



0(x, x) 



6 ^- 1-- 



x 



1 

1 - - 

z 



m c 



(54) 



where the derivative, dO{x,x)j dx, arises as a result of the integration by parts. Similar 
results are obtained for the contribution associated with the correlation function N(x,x), 
as well as for the corresponding contributions from (52) to be combined with 0(x, 0) and 
N(x,0). This indicates that all four nonperturbative functions of x, 0(x,x), N(x,x), 
0(x, 0) and N(x, 0), contribute both in the derivative and nonderivative forms to the twist- 
3 SSA, and the partonic hard cross sections convoluted with them are entirely determined 
from the hard scattering parts for the 7*g — >■ cc scattering which corresponds to the 2 — > 2 
process at the twist-2 level. 

Substituting (47) into (12) and using (50), (52)-(54), we obtain the leading-order QCD 
formula for the single-spin-dependent cross section Act in the SIDIS, ep^ — > eDX, generated 
from the twist-3 three-gluon correlation functions 0(xi,X2) and N(xi,X2) of (1) and (2) as 

d 6 Aa 



dxbj dQ 2 dzf dq\ d<\>d\ 



a 2 em a s e 2 c M N 
\^ 2 x%S 2 ev Q 2 



V Z 7 fc=l,-,4,8,9 X J Z 



1 - 



X 



1 - 



x sin($ 5 - x )Ak 



2q 7 



Q 2 



xa 



u 



dO(x,x) . 

x - 20(x,x) 

dx 



+ 



v z dq\ 



+ 



x 



m c 

l 2 ~Q 2 i 



dal 



z dx 



D(z) 



0(x, x) 



+ COs($ 5 - x) 

+ sin($ 5 - x) 



dA k a u k 



0(x,x) + (0(x,x) -> N(x,x)) 



d(f> zqT 

2q T \xa£ 



-2A k 



Q 2 



xx < do(x,o) ; 

x — -2O(x,0) 



z 



dx 
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+ 



'Q 2 da? x 
v z dq\ 



L/^ rl + 2A( _^ + ^- ) | o(xj0) 



zq T 



+ cos($ 5 - x) 



Q 2 



+ 



1-2 



1 — 2 <9x 



O(x,0) 



+ 



zq T 



-2 



+ (0(x, 0) ->■ -JV(x, 0)) 



(55) 



where a s = g 2 /(47r) is the strong coupling constant, Ak = Ak{<t> — x), and &l (j = 
U, XX, YY, {XY}) are defined in (53). It is worth noting the relations <7g 9 = for 
j = XX, YY and o"^ 2 34 = f° r 3 = {^^}? since even number of V M 's has to be in- 
volved in the contraction to give nonzero o\. We also remind the relations ^f- = — -4-8) 
= — 2.^9 = As and = 2A4 in (55). It is straightforward to derive the explicit 
formulae for a{ by calculating the diagrams for the 2 — > 2 processes in Fig. 3: For example, 



( = 



-XX 



cr^ Y ) are given in Eq. (81) of [19], where are nonzero for k = 1, • • • ,4, 
0. Those determine the gluon contribution to the twist-2 unpolarized cross 
section <j un P o1 for ep — > eDX, as 



while a$ g 



d 5 a unpo1 



dxbjdQ 2 dzfdqj>d<f) 87ixljS^ p Q 



y 2 \T,M<P)J d iJ d ii:D a (z)G( X )a 




(56) 



which can be obtained immediately from (22), (23) using (43) and (45), integrating over x 
for the fixed 0' = — x, and perfoming the formal replacement 0' — >■ (i.e., in (56) is 
understood to be this 0'). We do not show the rather lengthy formulae for the other hard 
cross sections in (53). 

Note that (55) shows the result for the c-quark fragmentation channel according to 
the diagrams in Fig. 3. The contribution for the c-quark fragmentation channel, due to 
the diagrams in Fig. 3 with the direction of the quark lines reversed, can be calculated 
similarly as above, and yields the formula (55) with the c-quark fragmentation function D(z) 
replaced by the c-quark fragmentation function and also with the replacements 0(x,x) — > 
—0(x,x) and 0(x, 0) — > — O(x,0), reflecting the fact that (1) is associated with the C-odd 
combination of the gluon operators. Combining this result with (55), we obtain the total 
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result, which can be recast into the following form: 
d 6 Aa 



dxbj dQ 2 dzf dq^dcpdx 
a 2 em a s e 2 c M N 



167r 2 4^g 2 V 2 ) k J^ AA9 AkSk l x I z S (q 2 i 1 x) i 1 z) + z 2 Q 2 j 



x E D a (z) 

a=c,c 



U d 20(x,x) 

MU 0( ^ } — x— 



)a^ + (A 0( ,,o)- 



x 



X 



f / d Ar/ , 2jV(x,x 



■)^i-(s^.o)- 



2AT(a;,0) > 



+ ivM A .3 _ jvm) A .; 

Ob Ob 



(57) 



where ^ = Ak{4> — X), and S k is defined as = sin($5 — x) f° r k — 1)2,3,4 and 
Sk = cos($5 — x) f° r fc = 8, 9. The quark-flavor index a can, in principle, be c and c, with 
<5 C = 1 and 5 5 = —1, so that the cross section for the Z)-meson production ep* — > eDX 
can be obtained by a simple replacement of the fragmentation function to that for the D 
meson, D a {z) — > D a (z). Comparing (57) with (55), we find 



Aal 



2qxx 



Q 2 {l-z) 



a 



k > 



(58) 



and, similarly, the partonic hard cross sections Aa k (j = 2,3,4, k = 1, • • -,4,8,9) are 
completely expressed by a k (j = U, XX, YY, {XY}) in (53). Substitution of the above- 
mentioned explicit forms of a k into these relations yields the formulae for Aa k ; these for- 
mulae, of course, agree with Eqs. (71)-(74) in [19], which were obtained by the direct 
calculation of the three-gluon diagrams in Fig. 2 in our previous study. 



4.3 Toward extension to higher orders 

As we have demonstrated in the last section, our master formula allows us to derive the 
explicit form of the twist-3 SSA for ep* — > eDX in the leading-order QCD, taking into ac- 
count the whole contribution induced by the three-gluon correlation inside the transversely- 
polarized nucleon. Thus, extension of our master formula beyond the leading-order QCD 
is interesting in that it would provide a powerful framework that allows us to calculate the 
higher-order corrections to the SSA, which are unknown at present. 
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To derive the master formula, our starting point was (13), which gives the whole twist-3 
cross section as the SGP contribution at x 1 = x 2 . The derivation of (13) presented in [19] 
relies on particular properties satisfied by the corresponding partonic scattering amplitudes 
with an on-shell internal line for which the propagator is replaced by its pole contribution. 
In particular, those properties include Ward identities for the relevant partonic amplitudes 
at the leading order. We expect that the similar Ward identities hold even after inclusion 
of the higher-order corrections employing a suitable gauge choice like the background-field 
gauge, and that (13) is useful beyond the leading order. Once the cross section is expressed 
as in (13), it is easy to see that the formulae (40), (42) hold for the SGP contributions: We 
note that the essential ingredient leading to these relations is the simple correspondence 
between the hard scattering part (31) relevant to the SGP contributions and the hard 
part (24) for the 7*g — > cc. As shown in (41) and (42), for the contributions that can be 
expressed as (26)-(31) using the 'jgcc vertices J 7 ^ and T\, the derivative with respect to fcj 
in (13) can be "transformed" into the derivative with respect to p\ in the master formula, 
without modifying J 7 ^ and T\. Similarly, for the higher-order contributions that can be 
expressed generically as (26)- (31) using the jg cc vertices J 7 ^ and T h p with the corresponding 
corrections included, the result (42) holds without additional modification to J 7 ^ and 

When including higher-order corrections, the SGP contributions may also occur from 
the diagrams that contain, in general, a number of final unobserved partons, as shown 
in Fig. 4. In that case, the vertex J 7 ^ in Fig. 4 has a number of external legs together 
with the loop correction inside. Still, we note that only the diagrams, which have the extra 
gluon-line attached to the final parton fragmenting into the D-meson, eventually contribute 
to the twist-3 cross section [2]. Because of this particular structure in the relevant SGP 
contributions, we can show that the relation (42) and thus our master formula hold for the 
twist-3 single-spin-dependent cross section for ep* — > eDX, even after the inclusion of the 
higher-order corrections. The details will be discussed elsewhere. 

5 Summary 

In this paper, we have derived the master formula for the contribution of the three-gluon 
correlation functions to the twist-3 single-spin-dependent cross section for the D-meson 
production in SIDIS, ep^ — > eDX. This formula connects the twist-3 effects due to the 
interference arising in the partonic hard cross section to the Born cross sections for the 
7*g — > cc scattering at the twist-2 level. In particular, the hard cross sections for the 
three-gluon correlation functions {0(x, x), N(x, x)} are completely determined by the hard 
cross sections associated with the gluon density distribution in the twist-2 unpolarized cross 
section for ep — > eDX. In the similar master formula derived for the contribution of the 
twist-3 quark-gluon correlation functions, only the SGP component of the cross section was 
connected to the twist-2 unpolarized cross section. For the present case with the three- 
gluon correlation functions, all contributions to the corresponding cross section appear as 
the SGP contribution, and thus the master formula derived here is for the total twist-3 cross 
section. The formula derived here can be easily extended to the three-gluon contribution 



21 



Figure 4: Generic diagrams for the higher-order corrections to the twist- 2 unpolarized 
cross section (upper diagram) and the twist-3 single-spin-dependent cross section (lower 
diagrams) . 

to p^p DX [25]. The derivation of the formula is based on the general structure and 
properties of the relevant twist-3 hard scattering part, which are expected to hold even 
after including the higher order corrections. 
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